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Abstract. We present an 0(n i log 2 rc)-time algorithm for the following problem: 
given a finite metric space X, create a star-topology network with the points of X 
as its leaves, such that the distances in the star are at least as large as in X, with 
minimum dilation. As part of our algorithm, we solve in the same time bound the 
parametric negative cycle detection problem: given a directed graph with edge 
weights that are increasing linear functions of a parameter X, find the smallest 
value of X such that the graph contains no negative-weight cycles. 



1 Introduction 

A metric space is a set of sites separated by symmetric positive distances that obey the 
triangle inequality. If X and Y are metric spaces and / :XhF does not decrease the 
distance between any two points, the dilation or stretch factor of / is 

d(f( Xl )J(x 2 )) ^ 
x u x 2 ex d(xi,x 2 ) 

We define a star metric to be a metric space in which there exists a hub h such that, for 
all x and y, d(x,y) — d(x,h) +d(h,y). Given the distance matrix of an n-point metric 
space X, we would like to construct a function / that maps X into a star metric Y, that 
does not decrease distances, and that has as small a dilation as possible. In this paper we 
describe an algorithm that finds the optimal / in time (9(n 3 log 2 n). Our problem may 
be seen as lying at the confluence of three major areas of algorithmic research: 

Spanner construction. A spanner for a metric space X is a graph G with the points of 
X as its vertices and weights (lengths) on its edges, such that path lengths in G equal 
or exceed those in X; the dilation of G is measured as above as the maximum ratio 
between path length and distance in X. The construction of sparse spanners with low 
dilation has been extensively studied [9] but most papers in this area limit themselves to 
bounding the dilation of the spanners they construct rather than constructing spanners 
of optimal dilation. Very few optimal spanner construction problems are known to be 
solvable in polynomial time; indeed, some are known to be NP-complete [15] and others 
NP-hard [3, 8]. Our problem can be viewed as constructing a spanner in the form of a 
star (a tree with one non-leaf node) that has optimal dilation. 

Metric embedding. There has been a large amount of work within the algorithms com- 
munity on metric embedding problems, in which an input metric space is to be em- 
bedded into a simpler target space with minimal distortion [16]; typical target spaces 



for results of this type include spaces with L p norms and convex combinations of tree 
metrics. As with spanners, there are few results of this type in which the minimum dila- 
tion embedding can be found efficiently; instead, research has concentrated on proving 
bounds for the achievable dilation. Our result provides an example of a simple class of 
metrics, the star metrics, for which optimal embeddings may be found efficiently. As 
with embeddings into low-dimensional L p spaces, our technique allows an input met- 
ric with a quadratic number of distance relationships to be represented approximately 
using only a linear amount of information. 

Facility location. In many applications one is given a collection of demand points in 
some space and must select one or more supply points that maximize some objective 
function. For instance, the 1 -median (minimize the sum of all distances from demand 
points to a single supply point) and 1 -center (minimize the greatest distance between 
any destination point and a single supply point) can be applied to operational challenges 
such as deciding where to build a radio transmitter or railroad hub so as to maximize 
its utility [7]. In a similar vein the problem discussed in this paper may be seen as 
selecting a single supply point to serve as the hub of a star-topology network. In this 
context dilation corresponds to the worst multiplicative cost penalty imposed on travel 
between any pair of input points due to the requirement that all travel is routed through 
the hub (center) point. Superficially, our problem differs somewhat from typical facility 
location problems in that the star we construct has a hub that is not given as part of the 
input. However, it is possible to show that the hub we find belongs to the tight span 
of the input metric space [6], a larger metric space that has properties similar to those 
of Loo spaces. Viewing our problem as one of selecting the optimal hub point from the 
tight span gives it the format of a facility location problem. 

Previously [10] we considered similar minimum dilation star problems in which 
the input and output were both confined to low-dimensional Euclidean spaces. As we 
showed, the minimum-dilation star with unrestricted hub location may be found in 
0{n\ogn) expected time in any bounded dimension, and for d = 2 the optimal hub 
among the input points may be selected in expected time 0(n2 a ^ log 2 n), where oc(n) 
is the inverse Ackermann function. For the general metric spaces considered here, the 
difficulty of the problems is reversed: it is trivial to select an input point as hub in time 
<9(« 3 ), while our results show that an arbitrary hub may be found in time (9(n 3 log 2 n). 

As we discuss in Section 2, the minimum dilation star problem can be represented 
as a linear program; however solving this program directly would give a running time 
that is a relatively high order polynomial in n and in the number of bits of precision of 
the input matrix. In this paper we seek a faster, purely combinatorial algorithm whose 
running time is strongly polynomial in n. Our approach is to first calculate the dilation 
X* of the optimal star. We do this by forming a X-graph G(X): a directed graph with 
weights in the form w(e) = X ■ m e + b e for parameters m e > and b e determined from 
the input metric. G(X) has the property that it contains no negative weight cycles if and 
only if there exists a star with dilation X. Next we calculate X*, the smallest value such 
that G(X*) contains no negative-weight cycles, which is also the dilation of the star we 
will eventually create. Finally we use G(X) and X* to compute the lengths of the edges 
from the star's center to each site, and output the resulting star. 



Our algorithm for computing X*, the smallest parameter value admitting no nega- 
tive cycles in a parametrically weighted graph, warrants independent discussion. To our 
knowledge no known strongly polynomial algorithm solves this problem in full gener- 
ality. Karp and Orlin [14] gave an 0(mn) time algorithm for a problem in which the 
edge weights have the same form w(e) = X ■ m e + b e as ours, but where each m e is re- 
stricted to the set {0, 1 }. If all m e = 1, the problem is equivalent to finding the minimum 
mean cycle in a directed graph [13], for which several algorithms run in 0(mn) time [4]. 
In our problem, each m e may be any nonnegative real number; it is not apparent how 
to adapt the algorithm of Karp and Orlin to our problem. Gusfield provided an upper 
bound [12] on the number of breakpoints of the function describing the shortest path 
length between two nodes in a A.-graph, and Carstensen provided a lower bound [2] for 
the same quantity; both bounds have the form n®( lo S"). Hence any algorithm that con- 
structs a piecewise linear function that fully describes path lengths for the entire range 
of X values takes at least n®( lo S") time. In Section 4 we describe our algorithm, which is 
based on a dynamic programming solution to the all pairs shortest paths problem. Our 
algorithm maintains a compact piecewise linear function representing the shortest path 
length for each pair of vertices over a limited range of X values, and iteratively contracts 
the range until a unique value X* can be calculated. Thus it avoids Carstensen's lower 
bound by finding only the optimal X*, and not the other breakpoints of the path length 
function, allowing it to run in (9(n 3 log 2 «) time. 



Fig. 1. Example of a metric space and its optimal star, which has dilation X* = 8/5. 



2 Linear Programming Formulation 

In this section we formally define the overall minimum dilation star problem and de- 
scribe how to solve it directly using linear programming. Our eventual algorithm never 
solves nor even constructs this linear program directly; however stating the underlying 
linear program and its related terminology will aid our later exposition. 

The input to our algorithm is a finite metric space. Formally, a metric space X is a 
tuple X = (X, dx) , where X is a set of sites and the function dx maps any pair of sites to 
the nonnegative, real distance between them. The following metric conditions also hold 
for anyx,y,z eX: 

1. dx(x,y) = if and only if x = y (positivity); 

2. dx(x,y) = dx(y,x) (symmetry); and 

3. dx(x,y) +dx{y,z) > dx{x,z) (the triangle inequality). 
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The input to our algorithm is a finite metric space S = (S,ds); we assume that the 
distance ds(x,y) between any x,y G S may be reported in constant time, for instance by 
a lookup matrix. 

A star is a connected graph with one center vertex. A star contains an edge between 
the center and every other vertex, but no other edges. Hence any star is a tree of depth 
1, and every vertex except the center is a leaf. Our algorithm must output a weighted 
star H whose leaves are the elements S from the input. The edge weights in H must 
be at least as large as the distances in S, and must obey reflexivity and the triangle 
inequality. In other words, if du(x,y) is the length of a shortest path from x to y in H, 
then d H {x,y) > ds{x,y), d H {x,y) = d H (y,x), and d H {x,y) + d H {y,z) > d H {x,z) for any 
vertices x,y,z in H. 

We also ensure that the dilation of H is minimized. For any two vertices u, v in some 
weighted graph G whose vertices are points in a metric space, the dilation between u 
and v is 

d G (u,v) 

8 G (",V) = —, r. 

ds(u,v) 

The dilation of the entire graph G is the largest dilation between any two vertices, i.e. 

Ag = max 8g(m, v). 

u.yeG 

Our output graph H is a star; hence every path between two leaves has two edges, so if 
we apply the definition of dilation to H, we obtain 

d H (u,c)+d H (c,v) w uc + w cv 
d H {u,v) = - 



d s (u,v) d s (u,v) 

where w XJ - is the weight of the edge connecting x and y in H. Hence the dilation of H 
may be computed by 

. W u . c + W c , v 

Ah = max ■ 



u,veH ds{u,c) 

This equation lays the foundation for our formulation of the minimum dilation star 
problem as a linear program. 

Definition 1. Let L be the following linear program, defined over the variables X and 
c v for every v <E S: 



such that for any v £ S, 
and for any v, w G S, 



Minimize X 

c v >0, (1) 

c v + c w > ds(v,w) (2) 

c v + c w < X ■ ds(v, w). (3) 



Let X* be the value assigned to X in the optimal solution to L. In other words, X* is the 
smallest dilation admitted by any set of distances satisfying all the constraints of L. 



L is clearly feasible. For example, if D = max Xyy£ sds(x,y), then the solution Vv c v =D 
and X = 2D/ mm Xtye sds(x,y) is a feasible, though poor, solution. 

Lemma 1. For any optimal solution of L, the value ofX gives the minimum dilation of 
any star network spanning S, and the c v values give the edge lengths of an optimal star 
network spanning S. 

Proof. Each variable c v corresponds to the weight w ViC of the edge between c and v in 
H. Inequality 1 ensures that the distances are nonnegative, Inequality 2 ensures that they 
obey the triangle inequality, and Inequality 3 dictates that A- is a largest dilation among 
any pair of sites from S. The value of X is optimal since L is defined to minimize X. 

Unfortunately L contains 0(n) variables and 0(n 2 ) constraints. Such a program 
could be solved using general purpose techniques in a number of steps that is a high- 
order polynomial in n and the number of bits of precision used, but our objective is 
to obtain a fast algorithm whose running time is strongly polynomial in n. Megiddo 
showed [19] that linear programs with at most two variables per inequality may be 
solved in strongly polynomial time; however our type (3) inequalities have three vari- 
ables, so those results cannot be applied to our problem. 

3 Reduction to Parameteric Negative Weight Cycle Detection 

In this section we describe a subroutine that maps the set of sites S to a directed, 
parametrically-weighted A.-graph G(X). Every edge of G(X) is weighted according to 
a nondecreasing linear function of a single graph-global variable X. An important prop- 
erty of G(X) is that the set of values of X that cause G(X) to contain a negative weight 
cycle is identical to the set of values of X that cause the linear program L to be infeasi- 
ble. Thus any assignment of X for which G(X) contains no negative weight cycles may 
be used in a feasible solution to L . 

Definition 2. A X-graph is a connected, weighted, directed graph, where the weight 
w(e) of any edge e is defined by a linear function in the form 

w(e) = X-m e + b e , 

where m e and b e are real numbers and m e > 0. 

Definition 3. Let G(X) be the X-graph corresponding to a particular set of input sites 
S. G(X) has vertices s and sfor each s G S. For s,t G S, G(X) has an edge of length 
—ds(s,t) from stot, andfors^=t, G(X) has an edge of length X-ds(s,t) from s to t_. 

Note that an edge from s to t has weight —ds(s,s) =0 when s = t . An example A.-graph 
G(X) for n = 3 is shown in Figure 2. 

Lemma 2. G(X) may be constructed in 0(n 2 ) time. 

Proof. G(X) has 2n vertices and 0(n 2 ) edges, each of which may be initialized in con- 
stant time. 



Fig. 2. The graph G(k) for n = 3. The weights of grayed edges are omitted. 



Lemma 3. IfX > 1 is assigned such that L has a feasible solution, then G(X) contains 
no negative weight cycle. 

Proof. Since G(X) is bipartite, any sequence of edges M traversed by a cycle in G(X) 
has even length. Depending on which partition M begins with, the sequence either takes 
the form 

M=((s-,s^),(s^,s-),(s-,sj ± ),...,(sj L1 s-)) 

or 

M = ((ijj , Si 2 ), [Si 2 , Si 3 ), (si 3 ,Sj 4 ), . . . , (si k ,S,'j )) , 

where i (1 ,i,- 2 , . . . ,si t are vertices from G(X). In either case, the cycle has weight 

w(M) = X ■ d s (s h , s h ) - rf s (s,- 2 , j i3 ) + X ■ d s (s i3 ,s k )-...- d s (s ik , s h ) (4) 

by the commutativity of addition. Since L is feasible, there exists some set of distances 
C satisfying the constraints of L, i.e. 

c x + c y <X- d s (x,y) (c x + c y )/X < d s (x,y) (5) 

and 

Cx + Cy > d s (x,y) => -(c x + c y ) < -d s (x,y). (6) 
Substituting (5) and (6) into (4), we obtain 

w(M) > X((c n + c h )/X) - (c h + c h ) + X((c, 3 +c k ))-...- (c ik + c n ) 

> (c h + c h ) - (c h + c h ) + (c, 3 + c /4 ) - . . . - (c k +c h ) 

^1 c ![ — + c !2 — C I2 + • • • + c i k — Ci k 

> 0. 

Theorem 1. Any set S ofn sites from a metric space may be mapped to a X-graph G(X) 
with 0(n) vertices, such that for any X > 1, G(X) contains a negative weight cycle if 
and only if L is infeasible for that value of X. The mapping may be accomplished in 
0(n 2 ) time. 

Proof. By Lemma 2, G(X) may be created in 0(n 2 ) time, and by Lemma 3, feasibility 
of L implies an absence of negative cycles in G(X). Section 5 describes an algorithm 
that, given a value X for which G(X) has no negative cycle, generates an edge length 
c v for every v e S that obeys the constraints of L. Thus, by the correctness of that 
algorithm, an absence of negative cycles in G(X) implies feasibility of L . 



4 Searching for X* 



We now turn to the problem of computing the quantity X*. This problem is an example 
of parametric negative weight cycle detection: given a A--graph G(X), find X* , the small- 
est value such that G(X*) contains no cycles of negative weight. Our algorithm func- 
tions by maintaining a range [X\, X2] which is known to contain X*. Initially the range 
is [— 00,+ooJ; over <9(logn) iterations, the range is narrowed until it is small enough 
that X* may be calculated easily. This approach is similar in spirit to Megiddo's general 
parametric search framework [17, 18], which, in loose terms, searches for the solution 
to an optimization problem by simulating the execution of a parallel algorithm for the 
corresponding decision problem. 

Our algorithm is presented in Listing 1. It is an adaptation of a parallel all pairs 
shortest paths algorithm based on matrix squaring [20]. The original algorithm uses a 
matrix D,(m,v), which stores the weight of the shortest path from u to v among paths 
with at most 2 ! edges. Each D, (m,v) may be defined as the smallest sum of two cells 
of Dj_i, and D[i g 2 n] defines the shortest paths in the graph. In the context of that orig- 
inal algorithm, edges and paths had real-number lengths, so it was sufficient to store 
real numbers in D,-. In the context of this paper, an edge's weight is a linear function 
of a variable X; hence the weight of a path is a linear function of X. Unfortunately the 
minimum-cost path between u and v may be different for varying values of X, so the 
weight of the shortest path from u to v is defined by the minima of one or more linear 
functions of X. Such a lower envelope of linear functions may be represented by a piece- 
wise linear function; hence each element of D, must store a piecewise linear function. 
Without further attention the number of breakpoints in these piecewise linear functions 
would grow at every iteration, and eventually operating on them would dominate our 
algorithm's running time. To address this, at every iteration we choose a new interval 
[X\ , X2] that contains no breakpoints, so that every D, may be compacted down to a 
single linear function. 

Lemma 4. For any X G [X\ , X2], the function Di(u, v) as computed in the listing evalu- 
ates to the weight of the shortest path from u to v among paths with at most 2' edges, or 
+°° if no such path exists. 

Proof. We argue by induction on ;'. In the base case i = 0, D,-(k,v) must represent the 
weight of shortest path from u to v that includes up to 2° = 1 edges. The only such paths 
are trivial paths, for which u = v and Dj(u, v) = 0, and single edge paths, for which the 
path length equals the edge length. 

For i > 1, each Dj(u, v) is first defined as the lower envelope of two entries of Z),_i in 
line 10, then redefined as a strictly linear function over the new smaller range [X\ , X2] in 
line 16, so we argue that the lemma holds after each assignment. In the first assignment, 
Dj(u, v) is defined to be the lower envelope of [A'-i ( u , w) + A'-i (w, v)] for all w € V; in 
other words, every w G V is considered as a potential "layover" vertex, and D;(w,v) is 
defined as a piecewise linear function that may be defined by differing layover vertices 
throughout the range , A.2] ■ By the inductive hypothesis, the values represent 
weights of minimum cost paths with at most 2' -1 edges; hence the resulting D, values 
represent weights of minimum cost paths with at most 2' -1 + 2 ! ~' = 2 ! edges. 



Listing 1 Computing the quantity X*. 
1: INPUT: A X-graph G(X) with n vertices V. 

2: OUTPUT: X*, the smallest value of X such that G(X) has no negative-weight cycles. 

3: Let X\ = — °° and Xi = 

4: INVARIANT: < 1* < X 2 

5: INVARIANT: D,(«,v) contains a linear function that represents the length of the shortest 
path from u to v among the subset of paths that use at most 2' edges, as a function of X, for 
any X € [Xi , X 2 ] 

6: Let Dq be an n x n matrix of piecewise linear functions. 

{0 if u = v 

X-m e +b e if G(X) contains an edge e from utov 
+co otherwise 
8: for i=l,2,..., [~log 2 n\ do 
9: for u, v e V do 

10: A(«,v) =min we v [£>,-_! (k, w) +£>,-_! (w,v)] 
11: end for 

12: Let B be the set of breakpoints of the piecewise linear functions stored in the entries of £),-. 

13: Perform a binary search among the values in B, seeking an interval bounded by two con- 
secutive breakpoints that contains X* . At each step, the test value of the binary search is 
less than X* if and only if setting X equal to the test value causes G(X) to contain a nega- 
tive cycle; use the Bellman-Ford shortest paths algorithm to determine whether this is the 
case. 

14: Set Xi and X 2 to the endpoints of the interval found in the previous step. 
15: for ii, v G V do 

16: Replace the piecewise linear function Di(u,v) with the equivalent linear function over 

the range [X1.X2]. 
17: end for 
18: end for 

19: Compute X*, the smallest value in the range [^A^], such that Z)j-(v,v) > for every v e V. 
20: Return X* . 



When D;(w,v) is reassigned in line 16, the range endpoints Ai and A2 have been 
contracted such that no entry of D, contains breakpoints in the range . Hence 

any individual D,(m,v) has no breakpoints in that range, and is replaced by a simple 
linear function. This transformation preserves the condition that D, -(u, v) represents the 
weight of the shortest path from u to v for any A G [A,i,A-2]. 

Lemma 5. Given two values Ai and A2 such that Ai < A-2, it is possible to decide 
whether X* < Ai, A* > A2, or A* G [Ai,A,2], in <9(n 3 ) time. 

Proof. By Lemma 3, for any value A', if G(A') contains a negative cycle when A = A', 
then A' < A*. So we can determine the ordering of Ai,A2, and A* using the Bellman- 
Ford shortest paths algorithm [1, 11] to detect negative cycles, as follows. First run 
Bellman-Ford, substituting A = A2 to evaluate edge weights. If we find a negative cycle, 
then report that A* > A2. Otherwise run Bellman-Ford for A = A4; if we find a negative 
cycle, then A* must be in the range [Ai,A2]. If not, then A* < A;. This decision process 
invokes the Bellman-Ford algorithm once or twice, and hence takes <9(n 3 ) time. 

Lemma 6. The algorithm presented in Listing 1 runs in 0(n 3 log 2 n) time. 

Proof. Each D,-_i (m, v) is a linear function, so each [D,-_i (u, w) + A'-i (w, v)] is a linear 
function as well. D,(m,v) is defined as the lower envelope of n such linear functions, 
which may be computed in <9(«log«) time [5]. So each D,(m,v) may be computed is 
0(n\ogn) time, and all 0(n 2 ) iterations of the first inner for loop take 0(n 3 logn) total 
time. Each D,(m,v) represents the lower envelope of 0(n) lines, and hence has 0(n) 
breakpoints. So the entries of D, contain a total of (9(n 3 ) breakpoints, and they may all 
be collected and sorted into B in (3(n 3 log«) time. Once sorted, any duplicate elements 
may be removed from B in (9(|-B|) = 0(n 3 ) time. 

Next our algorithm searches for a new, smaller [Ai , A2] range that contains A* . Recall 
that A* is the value of A for which G(A*) contains no negative weight cycle, and every 
entry of Z), is a piecewise linear function comprised of non-decreasing linear segments; 
so it is sufficient to search for the segment that intersects the A = line. We find this seg- 
ment using a binary search in B. At every step in the search, we decide which direction 
to seek using the decision process described in Lemma 5. Each decision takes (9(« 3 ) 
time, and a binary search through the 0(n 2 ) elements of B makes <9(logn) decisions, 
so the entire binary search takes <9(« 3 logn) time. 

Replacing an entry of D, with a (non-piecewise) linear function may be done naively 
in 0(n) time by scanning the envelope for the piece that defines the function in the range 
[Ai,A2]. So the second inner for loop takes (9(n 3 ) total time, and the outer for loop takes 
a total of 0(n 3 log 2 «) time. 

The initialization before the outer for loop takes 0{n 2 ) time. The last step of the 
algorithm is to compute A* , the smallest value in the range [Ai , A2] such that ( v, v) > 
for every v G V . At this point each D,-(h, v) is a non-piecewise increasing linear function, 
so this may be done by examining each of the n linear functions D^v, v), solving for its 
A-intercept, and setting A* to be the largest intercept. This entire process takes 0(n 2 ) 
time, so the entire algorithm takes <9(n 3 log 2 n) time. 

Theorem 2. The algorithm presented in Listing 1 calculates A* in <9(n 3 log 2 n) time. 



5 Extracting the Edge Weights 



Once X* has been calculated, all that remains is to calculate the weight of every edge 
in the output star. Our approach is to create a new graph G', which is a copy of G(X) 
with the addition of a new source node s with an outgoing weight edge to every 
v (see Figure 3). We then compute the single source shortest paths of G' starting at 
s, and define each c v to be a function of the shortest path lengths to v and v. This 
process is a straightforward application of the Bellman-Ford algorithm, and hence takes 
0(n 3 ) time. The remainder of this section is dedicated to proving the correctness of this 
approach. 



A* ■ d s (*3,»i) 




d s (si,s 



Fig. 3. The graph G' for n = 3. The weights of grayed edges are omitted. 



Definition 4. Let G' be a copy of the graph G(X) described in Definition 3, with all 
edge weights evaluated to real numbers for X — X*, and the addition of a source vertex 
s with an outgoing O-weight edge to every v £ G '. Let P(v) be a shortest path from s to 
v for any vertex v £ G ', and let l(v) be the total weight of any such P(v). The operation 
P(v)L)w yields the path formed by appending the edge (v, w) to P(v). 

Definition 5. Define c v = ^Hffl. 

We now show that our choice of c v satisfies all three metric space properties. 
Lemma 7. Every c v satisfies c v > 0. 

Proof. For each vertex v £ G' there exists an edge from v to v with weight 0. 

Lemma 8. Every distinct c v and c w satisfy c v + c w > ds(v, w) . 

Proof. By the definition of shortest paths, we have 

l(w) < l(v)—ds(v,w) 
ds(v,w) < l(v) — l(w). 



and by symmetric arguments, 

ds(w,v) < l(w) — l(v). 
Adding these inequalities, we obtain 

d s (v,w)+d s (w,v) < l(v)-l(w)+l(w)-l(v) 

ds(v,w) < (c v ) + (c w ). 

Lemma 9. Every distinct c v and c w satisfy c v + c w < X ■ ds(y, w) . 

Proof. Observe that the path P(w) U v is a path to v with weight /(w) +A,-rfj(w, v), and 
that the path P(v) Uw is a path to w with weight l(v) + X-ds(v,w). By definition P(v) is 
a shortest path to v, and similarly P(w) is a shortest path to w, so we have 

Z(v) <l(w)+X-d s (v,w) 

and 

/(w) < /(v)+X-d s (v,w). 
Adding these inequalities, we obtain 

l{v) + l(w) < (l(w) + \-ds(w,v)) + (l(v) + \-d s (v,w)). 
By assumption ds(w,v) — ds(v,w), so 

Z(v)-/(v)+/(w)-/(w) <2A,-*(v,w) 
(c v ) + (c w ) < k-ds(v,w). 

Theorem 3. Given S and the corresponding G(k) and X*, a set C of edge lengths c v for 
each v G S, such that for every v E S 

c v >0 

and for every distinct v,w E S 

c v + c w > d s (v,w) 
c v + c w <X-ds{v,w) 

may be computed in (9(n 3 ) time. 

Theorem 3 establishes that for any X* there exists a set C of valid edge lengths. This 
completes the proof of Theorem 1 . 



6 Conclusion 

Finally we codify the main result of the paper as a theorem. 

Theorem 4. Given a set S CI ofn sites from a metric space X — (X,d), it is possible 
to generate a weighted star H such that the distances between vertices of H obey the 
triangle inequality, and such that H has the smallest possible dilation among any such 
star, in 0(n 3 log 2 n) time. 



Acknowledgements This work was supported in part by NSF grant 0830403 and by 
the Office of Naval Research under grant N00014-08-1-1015. 



References 

1. R. Bellman. On a routing problem. Quarterly in Applied Mathematics, 16(l):87-90, 1958. 

2. P. J. Carstensen. Parametric cost shortest chain problem. Manuscript, Bellcore, 1984. 

3. O. Cheong, H. Haverkort, and M. Lee. Computing a minimum-dilation spanning tree is 
NP-hard. In Proc. 13th Australasian Symp. Theory of Computing, pages 15-24, 2007. 

4. A. Dasdan, S. Irani, and R. Gupta. An experimental study of minimum mean cycle 
algorithms. Technical report, University of California, Irvine, 1998. 

5. M. de Berg, M. van Kreveld, M. Overmars, and O. Schwarzkopf. Computational Geometry: 
Algorithms and Applications . Springer- Verlag, 2nd edition, 1998. 

6. A. W. M. Dress, K. T. Huber, and V. Moulton. Metric spaces in pure and applied 
mathematics. In Proc. Quadratic Forms LSU, Documenta Mathematica, pages 121-139, 
2001. 

7. Z. Drezner and H. Hamacher. Facility Location: Applications and Theory. Springer, 2002. 

8. J. Edmonds. Embedding into is Easy, Embedding into is NP-Complete. Discrete and 
Computational Geometry, 39(4):747-765, June 2008. 

9. D. Eppstein. Spanning trees and spanners. In J.-R. Sack and J. Urrutia, editors, Handbook 
of Computational Geometry, chapter 9, pages 425^-61. Elsevier, 2000. 

10. D. Eppstein and K. A. Wortman. Minimum dilation stars. Computational Geometry, 
37(1):27 - 37, 2007. Special Issue on the Twenty-First Annual Symposium on 
Computational Geometry - SoCG 2005. 

11. L. R. Ford and D. Fulkerson. Flows in Networks. Princeton University Press, 1962. 

12. D. Gusfield. Sensitivity analysis for combinatorial optimization. Technical Report 
UCB/ERL M80/22, Electronics Research Laboratory, Berkeley, 1980. 

13. R. M. Karp. A characterization of the minimum cycle mean in a digraph. Discrete 
Mathematics, 23(3):309-311, 1978. 

14. R. M. Karp and J. B. Orlin. Parametric Shortest Path Algorithms with an Application to 
Cyclic Staffing. Technical Report OR 103-80, MIT Operations Research Center, 1980. 

15. R. Klein and M. Kutz. Computing geometric minimum-dilation graphs Is NP-hard. In Proc. 
14th Int. Symp. Graph Drawing, volume 4372 of Lecture Notes in Computer Science, pages 
196-207, 2007. 

16. N. Linial. Finite metric spaces-combinatorics, geometry and algorithms. In Proc. 
International Congress of Mathematicians, Beijing, volume 3, pages 573-586, 2002. 

17. N. Megiddo. Combinatorial optimization with rational objective functions. In Proc. 10th 
ACM Symp. Theory of computing, pages 1-12. ACM, 1978. 

18. N. Megiddo. Applying Parallel Computation Algorithms in the Design of Serial 
Algorithms. / ACM, 30(4):852-865, 1983. 

19. N. Megiddo. Towards a Genuinely Polynomial Algorithm for Linear Programming. S1AM 
Journal on Computing, 12(2):347-353, 1983. 

20. C. Savage. Parallel Algorithms for Graph Theoretic Problems. PhD thesis, University of 
Illinois, Urbana-Champaign, 1977. 



